Probability distributions of Linear Statistics in Chaotic Cavities and associated phase transitions. 
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We establish large deviation formulas for linear statistics on the N transmission eigenvalues {Ti} of a 
chaotic cavity, in the framework of Random Matrix Theory. Given any linear statistics of interest A = 
X^i==i a (^i)> ' ne probability distribution 7a(A,N) of A generically satisfies the large deviation formula 
limjv->oo [— 2 log 7a{Nx, N)/f3N 2 ] — ^a(x), where ^>a(x) is a rate function that we compute explic- 
,-0 , itly in many cases (conductance, shot noise, moments) and /3 corresponds to different symmetry classes. Using 

^ ■ these large deviation expressions, it is possible to recover easily known results and to produce new formulas, 

such as a closed form expression for v(n) = limjv->oo var(T„) (where 7 n — 2~2i T?) for arbitrary integer 
£ — , n. The universal limit v* — lmin-joo v(n) = l/2n/3 is also computed exactly. The distributions display a 

central Gaussian region flanked on both sides by non-Gaussian tails. At the junction of the two regimes, weakly 
non-analytical points appear, a direct consequence of phase transitions in an associated Coulomb gas problem. 
Numerical checks are also provided, which are in full agreement with our asymptotic results in both real and 
Laplace space even for moderately small N. Part of the results have been announced in [P. Vivo, S.N. Majumdar 
and O. Bohigas, Phys. Rev. Lett. 101, 216809 (2008)]. 

C/3 ' 

. PACS numbers: 73.23.-b, 02.10. Yn, 21.10.Ft, 24.60.-k 

Keywords: Chaotic cavities, quantum dots, large deviations, Coulomb gas method 



C3 



O 

o 



- 1—^ 

X 



I. INTRODUCTION 



We consider the statistics of quantum transport through a chaotic cavity with Ni = N2 — N 1 open channels in the two 
attached leads. It is well established that the electrical current flowing through a cavity of sub-micron dimensions presents time- 
dependent fluctuations which persist down to zero temperature [1] and are thus associated with the granularity of the electron 
charge e. Among the characteristic features observed in experiments, we can mention weak localization O], universality in 
conductance fluctuations [3] and constant Fano factor fl. In the Landauer-Biittiker scattering approach 0], S la], the wave 
function coefficients of the incoming and outgoing electrons are related through the unitary scattering matrix S (2N x 2N): 
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where the transmission (i, t') and reflection (r, r') blocks are [N x N) matrices encoding the transmission and reflection coef- 
ficients among different channels. Many quantities of interest for the experiments can be extracted from the eigenvalues of the 
hermitian matrix ttf: for example, the dimensionless conductance and the shot noise are given respectively by G = Tr(#t) Jl 
andP = Tr[ttt(i_ tt t)] j^g. 

Random Matrix Theory (RMT) has been very successful in describing the statistics of universal fluctuations in such systems, 
and complements the fruitful semiclassical approach 1I9I- H3I1 . The simplest way to model the scattering matrix S for the case 
of chaotic dynamics is to assume that it is drawn from a suitable ensemble of random matrices, with the overall constraint of 
unitarity Il44l7l1 . Through a maximum entropy approach with the assumption of ballistic point contacts QJJ] , one derives that the 
prob ability distribution of S should be uniform within the unitary group, i.e. S belongs to one of Dyson's Circular Ensembles 

MM. 

It is then a non trivial task to extract from this information the joint probability density of the transmission eigenvalues {Ti} 
of the matrix W, from which the statistics of interesting experimental quantities could be in principle derived. Fortunately, this 
can be done 111 [l6ll20ll2lll and the expression reads: 



N 



P(T U ..., T N ) = A N J[ \Tj - T k f Ht^ 2 - 1 (2) 

j<k i—l 



where the Dyson index (3 characterizes different symmetry classes (/3 = 1,2 according to the presence or absence of time- 
reversal symmetry and /3 = 4 in case of spin-flip symmetry). The eigenvalues Tj are correlated random variables between and 
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1. The constant An is explicitly known from the celebrated Selberg's integral as: 



i JV TT ir ( 1 + i+^) r (f^ + 1 )) r ( 1 +^) 

A * = II b: . . ' r~77\ ^ 



fJo r(i + f)r(f + i + (iv + i-i)f) 



From ©, in principle the statistics of all observables of interest can be calculated. In this paper, we shall focus on the 
following: 



JV 



G = T l (conductance) (4) 

i=i 

N 

P = Y J T i( 1 - T i) ( shot noise ) (5) 

N 

T„ = T i (integer moments) (6) 



although in principle any linear statistics 1 A = Yli=i a (Ti) can be tackled with the method described below. It is worth 
mentioning that an increasing interest for the moments and cumulants of the transmission eigenvalues can be observed in recent 
theoretical |22] and experimental works |23Tl . 

Many results are known for the average and the variance of the above quantities, both for large N [Q],[24],[25l] and, very recently, 
also for a fixed and finite number of channels Ni, N2 B25I428I1 . In particular, a general formula for the variance of any linear 
statistics A = Yli=i a i^i) (where Ti = (1 + Ai) _1 ) in the limit of large number of open channels is known from Beenakker 
l29tl . However, at least for the case A — T„ (integer moments), this formula is of little practical use. The method we introduce 
below allows to obtain the sought quantity v(n) = lirriAr_ i . 00 var(T„) in a neat and explicit way. 

In contrast with the case of mean and variance, for which a wealth of results are available (see e.g. fl H^l and references 
therein), much less is known for the full distribution of the quantities above: for the conductance, an explicit expression was 
obtained for Ni = N2 = lj 2 JUrHU, while more results are available in the case of quasi one-dimensional wires l34ll and 3D 
insulators 05I1 . For the shot noise, the distribution was known only for Ni — N2 = 1 113611 . Very recently, Sommers et al. fl37ll 
announced two formulas for the distribution of the conductance and the shot noise, valid at arbitrary number of open channels 
and for any j3, which are based on Fourier expansions. Such results are then incorporated and expanded in l38ll . In l39ll and j40ll . 
along with recursion formulas for the efficient computation of conductance and shot noise cumulants, an asymptotic analysis 
for the distribution functions of these quantities in the limit of many open channels was reported. In a recent letter 14 ill , we 
announced the computation of the same asymptotics (in the form of large deviation expressions) u sing a Coulomb gas method 
and we pointed out a significant discrepancy with respect to the claims by Osipov and Kanzieper I39l l40ll . We will discuss in 
detail such disagreement and the way to sort it out convincingly in subsection llV Dl for the conductance case for f3 — 2. 

Given some recent experimental progresses ll42ll . which made eventually possible to explore the full distribution for the 
conductance (and not just its mean and variance), it is of great interest to deepen our knowledge about distributions of other 
quantities, whose experimental test may be soon within reach. 

It is the purpose of the present paper to build on Il4lll and establish exact large deviation formulas for the distribution of any 
linear statistics on the transmission eigenvalues of a symmetric cavity with N\ = N2 = N ^> 1 open channels. More precisely, 
for any linear statistics A whose probability density function (pdf) is denoted as 7a(A, N), we have 2 : 



lim 
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2\ogy A (A = Nx,N) 



PN 2 



= *a(x) (7) 



where the large deviation function ^a(x) (usually called rate function) is computed exactly for conductance (A = G), shot 
noise (^4 = P) and integer moments (A = T„). The method can be extended to the case of asymmetric cavities Ni = kN?, and 
is based on a combination of the standard Coulomb gas analogy by Dyson and functional methods recently exploited in B43I1 in 
the context of Gaussian random matrices. This approach has been then fruitfully applied to many other problems in statistical 
physics HQ. 



1 A is a lineal' statistics in that it does not contain products of different eigenvalues. The function a(x) may well depend non-linearly on x. 

2 Hereafter we will use the concise notation T(A, N) ^ exp ( — ^N 2 ty a{A/N) J to mean exactly (7). 
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The paper is organized as follows. Section II provides a quick summary of our main results. In Section III, we summarize 
the Coulomb gas method that can in principle be used to obtain the rate function associated with arbitrary linear statistics. In 
Section IV, V and VI, we use this general method to obtain explicit results respectively for the conductance, shot noise and 
integer moments. Finally we conclude in Section VII with a summary and some open problems. 



II. SUMMARY OF RESULTS 



Distribution of conductance G: We obtain the following exact expression for the rate function in the case of conductance (G) 

{I - log(4x) for < x < \ 

8(x-i) 2 for \<x<\ (8) 

| - log[4(l - x)] for | < x < 1 

The rate function has a quadratic form near its minimum at x = 1/2 in the range l/4<£<3/4. Using it then follows that 
the distribution ?q(G, N) has a Gaussian form close to its peak 



7 G {G,N) «exp 



--N z -8 

N 2 



G 1 



exp 



1 ( c N 
2(1/8/3) ^ 2 



(9) 



from which one easily reads off the mean and variance (G) = N/2 and var(G) = 1/8/3, which agree with their known 
large N values. The fact that the variance becomes indepedent of N for large N is referred to as the universal conductance 
fluctuations. This central Gaussian regime is valid over the region N/4 < G < 3N/4. Outside this central zone, J > G (G, N) has 
non-Gaussian large deviation power law tails. Using $Q in © near x = and x = 1 , we get 9 G (G, N) ~ G pN I 2 (as G ->• 0) 
and J , g(G, N) ~ (N — G)^ N / 2 (as G — > N) which are in agreement, to leading order in large N, with the exact far tails 
obtained in IL37L L38L l52h . It turns out that while the central Gaussian regime matches continuously with the two side regimes, 
there is a weak singularity at the two critical points G = N/4 and G = 3N/4 (only the 3rd derivative is discontinuous). One of 
our central results is to show that these two weak singularities in the conductance distribution arise due to two phase transitions 
in the associated Coulomb gas problem where the average charge density undergoes abrupt changes at these critical points. 
We note that an intermediate regime with an exponential tail claimed in l39l l4(ill does not appear in our solutions, and we will 
comment in detail about such discrepancy in subsection llV Dl 

Distribution of shot noise P: For the shot noise, the exact rate function reads 



±-21og2-±log 



* p(a:) = <64(x-±) 



for < x < 4s 

— — 16 

for Ta < X < A 
16 — — 16 



(10) 



i - 21og2 - £ log (i - x) for ^ < x < \ 



Thus the shot noise distribution CPp(P, N), in the large N limit, also has a central Gaussian regime over N/16 < P < 3N/16 
where 



y P (P,N) w exp 



B , ( P 



exp 



B „ (PI 



exp 



1 



2(1/64/3) 



P — 



(11) 



yielding (P) = N/8 and var(P) = 1/64/3. As in the case of conductance, this central Gaussian regime is flanked on both sides 
by two non-Gaussian power law tails with weak third order singularities at the transition points P = N/16 and P = 3N/16. 

Distribution of integer moments T„ : For general n > 1, the computation of the full rate function ^<y n (x) and hence the large 
N behavior of the distribution 3V„ N) turns out to be rather cumbersome. However, the distribution 3V n (^n, N) shares 
some common qualitative features with the conductance and the shot noise distributions. For example, we show that for any 
n > 1, there is a central Gaussian regime where the rate function $!<x n (x) has a quadratic form given exactly by: 



r(n+l/2) 



(12) 
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where b n = 4 [r(n+i%y\^ ■ This implies a Gaussian peak in the distribution: 



3V n (T„,A0«exp 



,b 

2 2 



T 



r(n + l/2) 



n 2 



AT 0Fr(n + l)_ 
from which one easily reads off the mean and variance of integer moments: 



(13) 
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FIG. 1: (Color online). ^g{x) and $p(s) as in ^ and JlOt . The black dots highlight the critical points on each curve. 



, s AT(n+l/2) 
y/nT(n + 1) 

To the best of our knowledge, the latter result, together with its universal asymptotic value linin^co var(T„) = (2(3tt)~ 1 , has 
not been reported previously in the literature. 

As in the case of conductance and shot noise, there are singular points separating the central Gaussian regime and the non- 
Gaussian tails. However, the situation is more complicated for n > 1. Unlike in the conductance or the shot noise case where 
there are only two phase transitions separating three regimes, there are more regimes for n > 1. We analyse here in detail the 
case n = 2 where we show that there are actually three singular points separating 4 regimes. This is again a consequence of 
phase transitions in the associated Coulomb gas problems where the average charge density undergoes abrupt changes at these 
critical points. It would be interesting to see how many such phase transitions occur for n > 2. However, we do not have any 
simple way to predict this and this remains an interesting open problem. The tails of the distribution for the integer moments 
case can be computed in principle for general n > 1, and the formalism is developed in section[VT] but for clarity we will focus 
mostly on the case n = 2. The left tail has again a power-law decay for any n > 1. 

III. SUMMARY OF THE COULOMB GAS METHOD 

It is useful to summarize briefly the method we use to compute the distribution of linear statistics in the large A^ limit. Given 
a linear statistics of interest A — a (Ti), its probability density function (pdf) "Pa(A, N), using the joint pdf of Tj's in Eq. 

©, is given by 

y A (A, N)=A N jf ■ • • jf dTi • • • dT N exp ^ lo § \ T i ~ T «\ + (f - X ) E lo S T *j 6 (j2 a ^ - A ) ■ < 16) 



The main idea then is to work with its Laplace transform 



poo 

F N (z;A) = / dAy A {A.N)e 
Jo 



- zA 



(17) 



where the variable z takes in principle complex values and its dependence on A^ is at present unspecified. 
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Suppose one is able to compute the following limit: 

t ( \ 2 r logFjv (N a p;A) 

J A (p) = -^hm W2 (18) 

and such limit is finite and nonzero for a certain speed a and for p G R such that p ~ 0(1) for iV — > oo. It is then a classical 
result in large deviation theory (Gartner-Ellis theorem, see e.g. [53] Appendix C) that the following finite nonzero limit exists: 

m f , 2 \ogV A (N a x, N) 

^) = -- jv lrm o W2 (19) 

and ^a(x) [the rate function] is given by the inverse Legendre transform of Ja(p)- 

^a(x) = maxp [Ja(p) —px]. (20) 

Note that: 

1 . There is no need to consider complex values for z (and thus for p) in ( fTTI ). as only real values matter for obtaining the rate 
function. As 7a(A, N) has generically a compact support, negative values forp are also allowed. 

2. Setting the appropriate speed a is evidently crucial for obtaining a finite nonzero limit in ( TT8l . There is no freedom in here. 
Conversely, any speed is equally good when extracting the cumulants out of the Laplace transform ( fTTI i (z = N a p) through 
the formula (set f3 — 2 for simplicity): 



1 (^r) ^lo g F N (Ny;A) 



(21) 

p=0 



as no limiting procedure is performed in (fJTJ. 



3. The rate function ^a(^) encodes the full information about the probability distribution in the limit of infinitely many open 
channels. However, our numerical simulations confirm that it also gives a fairly accurate description of such distribution for 
rather small N. 

What is then the correct speed a for the linear statistics considered here? It is quite easy to argue that a must be set equal to 1. 
The reason is best understood by taking the conductance G — YliLi Tt for /3 = 2 as an example. By very general arguments 
we expect the large N behavior of 7g(G, N) to scale as CPg(G, N) ~ e~ N *g(G/jv) f or j^gg _/y_ Clearly, the two exponentials 
(the one coming from ?g(G, N) and the other coming from the Laplace measure) must be of the same order in N to guarantee 
a meaningful saddle point contribution, and since G ~ N for large N, clearly z ~ N as well. 
After setting the proper speed, we get in full generality: 



V A (A, N)e~i NpA dA = A N dT x ■ ■ ■ dT N ex P ( f E lo S l T i " T *\ + (f " X ) S logTi " f^ZX 2 *) 



F N (Np;A) 



(22) 

We can write the exponential as, exp [-0E({Ti})] with E({Ti}) = -(1/2) V . , fc log |Tj - T fc | + £\ V(Tj) where V(T) = 
(1//3 — 1 /2) log(T) + pNa(T) /2. This representation provides a natural Coulomb gas interpretation. We can identify Tj's as 
the coordinates of the charges of a 2-d Coulomb gas that lives on a one dimensional real segment [0, 1]. The charges repel each 
other via the 2-d logarithmic Coulomb potential and in addition, they sit in an external potential V(T). Note that the Laplace 
parameter p appears explicitly in the external potential V(T). Then E is the energy of this Coulomb gas. Thus one can write the 
Laplace transform as the ratio of two partition functions 

F,v (Np; A) = ^ y A (A, N)e~^ NpA dA = §^3 (23) 
Jo AK-W) 



where Z p {N) is precisely the multiple integral on the rhs of Eq. d22l > and Zq(N) — 1/Ajy (which simply follows by putting 
p = in Eq. d22i > and using the fact that the pdf 7a(A, N) is normalized to unity). 

The next step is to evaluate this partition function Z p (N) of the Coulomb gas in the large N limit. This procedure for the 
large N calculation was originally introduced by Dyson 019n and has recently been used in the context of the largest eigenvalue 
distribution of Gaussian 14311 and Wishart random matrices [51] and also in other related problems of counting stationary points 
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in random Gaussian landscapes Il44tl . There are two basic steps involved. The first step is a coarse-graining procedure where 
one sums over (partial tracing) all microscopic configurations of Tj's compatible with a fixed charge density function g p (T) = 
N^ 1 J2i ${T — Ti) and the second step consists in performing a functional integral over all possible positive charge densities 
Qp(T) that are normalized to unity. Finally the functional integral is carried out in the large N limit by the saddle point method. 
Following this general procedure summarized in B43I1 . the resulting functional integral, to leading order in large N , becomes: 



Z p (N) oc j X>[ft,]e-** s[ ^ ] (24) 

where the action is given by 

S[Qp] = V [ Qp(T) a{T)dT + b\ f Q p (T)dT - 1 
Jo Uo 

- [ [ dTdT' g p {T)g p {T')\og\T -T'\. (25) 
Jo Jo 



Here B is a Lagrange multiplier enforcing the normalization of g p (T). In the large N limit, the functional integral in 
is particularly suitable to be evaluated by the saddle point method 3 , i.e., one finds the solution g p (T) (the equilibrium charge 
density that minimizes the action or the free energy) from the stationarity condition SS[g p ]/5g p = which leads to an integral 
equation 



V cxt (T) =pa(T) + B = 2 f Qp(T') log |T - T'\ dT' 

Jo 



(26) 



where V cx t(T) = pa(T) + B is termed as external potential. Differentiating once with respect to T leads to a singular integral 
equation 



^a'(T) = Pr C ^-IdT 1 
2 V ; J T-T' 



(27) 



where Pr denotes the principal part and a'(T) = da/dT. 

Assuming one can solve d27l ) for g*, one next evaluates the action S[g p ] in d25l ) at the stationary solution g* and then takes 
the ratio in ( 1231 to get (upon comparison with (fT8ll): 



3 A <P) 



V A (A,N)e~? N i> A dA » e~% N2[S ^ ~ (28) 



F N (Np;A) 



Inverting the Laplace transform gives the main asymptotic result "P(A, N) w cxp (~^N 2 ty a (^)J where the rate function is 
the inverse Legendre transform (see d20l>). 

^a(x) = m&x[-xp + Ja(p)] (29) 

p 



with Ja(p) given by the free energy difference as in I 

To summarize, given any linear statistics a(T), the steps are: (i) solve the singular integral equation (|27| | for the density g p {T) 
(ii) evaluate the action S[g*] in (l25l l (iii) evaluate Ja(p) = S[g*] — S[go\ and finally (iv) use Ja{p) in d29l , maximize the rhs 
to evaluate the rate function $ a(x). We will see later that all these steps can be carried out fully and explicitly when A = G 
(conductance) and A = P (shot noise) and partially when A — T„ (integer moments). 

The important first step is to find the explicit solution of the singular integral equation (|27| |. Note that this equation is of the 
Poisson form and it is, in some sense, an inverse electrostatic problem: given the potential pa'(T), we need to find the charge 



3 Note that such a nice feature is a direct consequence of having employed the correct speed a = 1, i.e. of having scaled the Laplace parameter z with N in 
the correct way. 
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density g*(T). To proceed, we recall a theorem due to Tricomi [54] concerning the general solution to singular integral equations 
of the form 



g(x) = Pr f — — —dx' 

J a ^ % 



(30) 



where g(x) is given and one needs to find f(x) which has only a single support over the interval [a, b] with a < b (the lower 
edge a should not be confused with the linear statistics function a(T)). The solution f(x), with a single support over [a, b] can 
be found explicitly B54I1 



/(*) 



n 2 y/(b — x){x — a) 



b y /(b-x')(x'-a) 
Pr / g(x )dx + B\ 

x ~ x' 



(31) 



where Bi is an arbitrary constant. 

In our case, g(x) = pa'(x)/2 and provided we assume that the charge density has a single support over [a, b] with a < b, 
we can in principle use this solution OTb . However, if the solution happens to have a disconnected support one cannot use this 
formula directly. Whether the solution has a single or disconnected support depends, of course, on the function a(T). We will 
see that indeed for the case of conductance (a(T) = T), the solution has a single support and one can use ( f3TT > directly. The 
edges a and b in that case are determined self-consistently as explained in Section IV. On the other hand, for the shot noise 
(a(T) = T(l — T)) and for integer moments with n > 1 (a(T) = T n ), it turns out that for certain values of the parameter p, 
the solution has a disconnected support. In that case, one cannot use ( f3TT > directly. However, we will see later that one can still 
obtain the solution explicitly by an indirect application of OTb . 

A very interesting feature of (f2Tb is that, depending on the value of the Laplace parameter p, the fluid of charged particles 
undergoes a series of real phase transitions in Laplace space, i.e., as one varies the Laplace parameter p, there are certain critical 
values of p at which the solution Qp~{T) abruptly changes its form. As a consequence, the rate function, related to the Laplace 
transform via the Legendre transform ( |29l l, also undergoes a change of behavior as one varies its argument at the corresponding 
critical points. The rate function is continuous at these critical points but it has weak non-analiticities (characterized by a 
discontinuous third derivative). 






FIG. 2: Phases of the density of transmission eigenvalues for the conductance case. 



As an example, we consider the case of the conductance (a(T) = T) (see section [IV] for details). In fig. |2j we plot schemat- 
ically the saddle point density Qp(T) (solution of ( 1271 ) with a(T) = T), for three different intervals on the real p line. We will 
see in the next section that there are three possible solutions valid respectively for p > 4, —4 < p < 4 and p < — 4. 

• when p > 4, the external potential V ex t (T) = p T + B is strong enough (compared to the logarithmic repulsion) to keep the 
fluid particles confined between the hard wall at T = and a point L\ = 4/ p. The gas particles accumulate towards T — > + , 
where the density develops an inverse square root divergence, while Qp(Li) = 0. This situation is depicted in the left panel of 
fig- El 
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• when p hits the critical value = 4 from above, the density profile changes abruptly. The external potential V cx t is 
no longer overcoming the logarithmic repulsion, so the fluid particles spread over the whole support (0,1), and the density 
generically exhibits an inverse square root divergence at both endpoints (T — > + and T — > 1~). This phase keeps holding 
for all the values of p down to the second critical point p(~> = —4 (see the second panel in fig. |2|, when the negative slope of 
the potential is so steep that the particles can no longer spread over the whole support (0, 1), but prefer to be located near the 
right hard edge at T = 1. 

• In the third phase (p < —4), the fluid particles are pushed away from the origin and accumulate towards the right hard wall at 
T — > 1 (see the rightmost panel in fig. [2]). The density thus vanishes below the point L 2 = 1 — 4/|p|. 

It is worth mentioning that such phase transitions in the solutions of integral equations have been observed recently in other 
systems that also allow similar Coulomb gas representations. These include bipartite quantum entanglement problem 15511 . non- 
intersecting fluctuating interfaces in presence of a substrate j47] and also multiple input multiple output (MIMO) channels l56ll . 

IV. DISTRIBUTION OF THE CONDUCTANCE 

We start with the simplest case of linear statistics, namely the conductance G = X)i=i ^i- Thus in this case a(T) = T is 
simply a linear function. Substituting a(T) = T in d26l l gives 

WT) = P T + B = 2 f g* p (T')log\T-T'\dT' (32) 
Jo 

and then ( |27| i becomes 

We have then to find the solution to P3t . Once this solution g*(T) is found, we can evaluate the action S[g*] at the saddle 
point in the following way. Multiplying ( 1321 by g*(T) and integrating (using the normalization J Q g*(T)dT = 1) gives 

pf T Q ;(T)dT + B = 2 f [ g;(T)g;(T r )x 
Jo Jo Jo 

x \og\T -T'\dTdT'. (34) 

Next we use this result to replace the double integral term in the action in d25l ) (with a(T) = T) to get 

S[q*] = ^J e ;(T)TdT-^. (35) 

The yet unknown constant B is determined from d32b upon using the explicit solution, once found. 

To find the solution to d33l ) explicitly we will use the general Tricomi formula in OTb assuming a single support over [a, b]. 
The edges a and b will be determined self-consistently. Physically, we can foresee three possible forms for the density Q*(T) 
according to the strength and sign of the external potential V ex t (T) = P T + B on the left hand side (lhs) of d32l . 

1 . For large and positive p, the fluid particles (transmission eigenvalues) will feel a strong confining potential which keeps them 
close to the left hard edge T -t 0+. Thus, g*(T) will have a support [0, L x ], with < L x < 1. 

2. For intermediate values of p, the particles will spread over the full range [0, 1]. 

3 . For large and negative p, the fluid particles will be pushed towards the right edge and the support of g* (T) will be over [L 2 , 1] , 
with < L 2 < 1. 

These three cases will correspond to different solutions for the Tricomi equation d33l above, and the positivity constraint for 
the obtained densities will fix the range of variability for p in each case. Once a solution gp~(T) for each case (different ranges 
for p) is obtained, we can then use ((35) to evaluate the action at the saddle point. 

Let us consider the three cases discussed above separately. 
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A. Large p: support on [0, Li] 

We assume that the solution is nonzero over the support [0, L\\ based on our physical intuition for large p, where L\ is yet 
unknown. We use the general Tricomi solution with a single support in ( f3Tb with a = 0, b = Li and g(x) — p/2 giving 



Q* P (T) ■ ' 



lfiy/T{L X -T) 



2 L T- T> 



P 



+Bj] (36) 
where B\ is an arbitrary constant. Evaluating the principal value integral on the rhs of (l36*l l we get 

^) = ^ 7T VT^T (37) 

where the constant B\ has been determined using the fact that the density gp~(T) must vanish at the upper edge T = L±. The 
normalization of Qp(T) gives: 

4 

L-l = -, Li<l^p>4 (38) 
P 

As expected, this solution holds for large values of p (i.e. for a strong confining potential). 

Since the point T = belongs to the support [0, Li], we can put T = in d32l to determine the constant B 

B = 2 f Q*(T') log T'dT' (39) 
Jo 

Substituting B in ( f35l > gives the saddle point action 

S[Q$] = | f Q* P (T)TdT - f e ;(T) \og(T)dT (40) 
1 Jo Jo 

Performing the integrals using the explicit solution q* (T) d37l i gives a very simple expression, valid for p > 4, 

S[g* p ] = 3/2 + logp (41) 

In Fig. |3]we show the results from a Montecarlo simulation to test the prediction (f3Tb for the average density of eigenvalues 
in Laplace space for p > 4. The numerical density for N eigenvalues (here and for all the subsequent cases) is obtained as: 

( e -^^-n i<fc K--^i 2 ) w , 

Q*M) * (42) 

(a-P»^*<Il j<k \xs-Xk\*) N 

where the average ^ • ^ is taken over N — 1 random numbers x%, . . . , Xn (numerator) with a flat measure over [0, 1], with 

xi spanning the interval (0, 1). In the denominator, the normalization constant is obtained with the same procedure, this time 
averaging over N random variables uniformly drawn from (0, 1). In all cases, the agreement with the theoretical results is fairly 
good already for N = 5. 



B. Intermediate p: support on the full range [0, 1] 

In this case, the solution of d33l from OTb reads: 

q*(T) = P [K - T] (43) 

pV ' 2vr v /T(l - T) 1 J 

The normalization of Qp(T) determines K = (4 + p) /2p. Now, depending on whether p > or p < 0, there are 2 positivity 
constraints {g*(T) > everywhere) to take into account: 

1. if p > 0, the positivity constraint K — 1 > at the upper edge T = 1 implies p < 4. 
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FIG. 3: (Color online). Density of transmission eigenvalues T for N = 5 and p = 6 (theory vs. numerics) for the conductance case. 



2. if p < 0, the positivity constraint K < at the lower edge T = implies p > —4. 

Thus the solution d43b with support over the full alllowed range [0, 1] is valid for all —4 < p < 4. 
Substituting this solution into the simplified action d40b (which holds in this case as well) gives: 

S[$ = ~ + f +21og2 (44) 

Note that since this range —4 < p < 4 includes, in particular the p = case, we can use the expression in ( |44T > to evaluate the 
value of the action at p = that will be required later in evaluating the large deviation function via (|281 l. Putting p = in (l44l > 

gives 



S[g* ] =2 log 2. 

In fig. |U we plot the analytical result for the density together with Montecarlo simulations for N = 5 and p = 1. 



(45) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 09 



FIG. 4: (Color online). Density of transmission eigenvalues T for N — 5 and p = 1 (theory vs. numerics) for the conductance case. 



C. Large negative p: support on [L2, 1] 




0.1 0. 



0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 5: (Color online). Density of transmission eigenvalues T for TV = 5 and p = —6 (theory vs. numerics) for the conductance case. 



In this case, the solution of d33t reads: 



\P\ 



\/T — L 2 , L 2 



4 

1_ ¥\ 



(46) 
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The implications for the p-range are as follows: 



L 2 > =>• p < -4 



(47) 



Note that we can no longer use the expression for the constant B in d39b since now the allowed range of the solution [La, 1] 
does not include the point T = 0. Instead, putting T = 1 in d32l . we determine the value of B as 



(48) 



B = 2 g*(x')log\l-x'\dx' -p 
Jo 

Substituting B in d35l > we get a new expression for the action at the saddle point 

S[q*] = | f dTT g*(T) - J 1 dT g ;(T) log |1 - T\ + | 

Evaluating ( |49l using the solution in (|46*T > gives the saddle point action for p < —4 

S[g* p ] -3/2+p + log(-p) 

In fig. [5] we plot the analytical result for the density together with Montecarlo simulations for TV — 5 and p = —6. 

logfF^Np, G)) 



(49) 



(50) 



w 2 




FIG. 6: (Color online). Left: Jg(p) and Jq(p) vs. Montecarlo simulations (see eq. d70b and d80ll). Right: our asymptotic predictions Jg(p) 
(solid red line) is compared with the exact finite N expression d59t from 13911 for N — 4 (green squares) and iV = 10 (blue dots). Already for 
TV = 4, our large TV formula again matches the exact finite TV result with accuracy up to the second decimal digit over the full real p range. 



D. Comparison with other theories and numerical simulations 

As summarized in the next section, we obtained for the rate function of the conductance for f3 — 2, defined as (see JT9]l): 

T , , logIP G (TVx,TV) 

* G (z) = - im -2 ^ >- (51) 

AT->oo TV 

the following expression (limited to x € [1/2, 1], given the symmetry ^g(x) = ^g(1 — x)): 

v 2 



1/2 < x < 3/4 
Dg[4(l - x)] 3/4 < sc < 1 



In [39], Osipov and Kanzieper (OK) claim a different limiting law, namely: 



*ok (x) = (8(*-D 2 l/2<*<3/4 



4x-| 3/4<2:<l 



(53) 
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and ^q K (x) would approach the form in (52\ (second line) only at the extreme edge x —} 1 (over a narrow region of order 
l/N). 

Which law is then correct? There is a conclusive way to settle this dispute, namely to compare the two theoretical results to a 
direct numerical simulation of 7c{G, TV). We will present simulation results in Laplace space which agree very well with our 
result on the Laplace transform (see fig. [6] left panel, and equations (l65l l and ( F70b in next section) over the full range of real p 
values, as well as a convincing comparison with the exact finite TV result for the same observable using the Hankel determinant 
representation (|59l from [39] (see fig. [6j right panel). We shall argue below that OK asymptotic theory is instead unable to 
reproduce the tails of Jg(p) for \p\ > 4, which are responsible for long power-law tails in the rate function ^g( x )- Since, 
however, working in the Laplace space may not appear conclusive as far as the real-space rate function ^g(x) is concerned, it 
would be better if one could perform a simulation directly for '5 > q(G 1 TV) and not just for its Laplace transform. 

Indeed, it turns out to be quite easy to simulate directly Vq(G, TV) using an elementary and standard Monte Carlo Metropolis 
algorithm which we describe below. 

Monte Carlo method: 

The main problem is to compute the distribution of the conductance G which, for a fixed number of channels TV and (3 — 2, 
is given by the multiple integral 

! N /N \ 

T G (G,N)=A N HdTiHiTi-TjfS ^Ti-G (54) 

•'O i—l i<j \i=l / 

where the prefactor An is set by the normalization: J °° ?q{G, N)dG = 1 and is known exactly for all TV (see eq. (01). 




G/N x=G/N 

FIG. 7: (Color online). Comparison between Montecarlo simulations for the rate function = — limjv^oo log Tg^x, N)/N 2 as a 

function of % = G/N for ft = 2 (black dots), our asymptotic theory (VMB) (continuous red line) and the theory by Osipov-Kanzieper (OK) 
(dashed blue line). Simulations are performed for TV = 10 (left) and TV = 50 (right), which allows to appreciate the convergence to our 
theoretical curve as TV is increased. 



To employ the Monte Carlo method, we first write the integrand (the Vandermonde term) in d54l ) 

nm-T,) 2 ^-^}; E{T i } = -^(\T i -T j \). (55) 

i<j ijtj 

Thus, one can interpret < Tj < 1 as the position of an i-th charge in a one dimensional box [0, 1] and the charges interact via 
the logarithmic Coloumb energy — E{Ti]. This Coulomb gas is in thermal equilibrium with a Gibbs weight cxp[— 7?{Tj}] for 
any configuration {T;}, where the inverse temperature is set to 1. 

It is then very easy and standard to simulate the equilibrium properties of this gas via a Monte Carlo method {5% . We start 
from any configuration {T^}. We pick a particle, say the i-th one, at random and attempt to move its position by an amount S: 
Ti — > Ti + 5. This move causes a change in energy AE of the gas. According to the standard Metropolis algorithm lf57ll . the 
move is accepted with probability e~ AE if AE > and with probability 1 if AE < 0. The move is rejected if the new position 
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Ti + S is outside the box [0, 1]. This Metropolis move guarantees that after a large number of microscopic moves, the system 
reaches the stationary distribution with the correct Boltzmann weight e~ E ^ Ti ^ . 

We wait for a long enough time to ensure that the system has indeed reached equilibrium. After that we let the system evolve 
according to these microscopic moves and construct the normalized histogram 7q(G, N) of G — YliLx ^i- Once again, we 
are guaranteed that TVs are sampled with the correct equilibrium weight. This procedure allows us to simulate fairly large 
systems. To obtain good statistics for the distribution over the full range of a; = G/N, i.e., over 1/2 < x < 1, we implement an 
iterative conditional sampling method, used in other contexts before [49], that allows us to generate events with extremely small 
probabilities at the far tail of the distribution I158ll . 

In Fig. H we plot $(x) = - ln[7 G (G, N)]/N 2 vs the scaling variable x = G/N for 1/2 < x < 1 for N = 50. The 
black dots show the simulation points. The red (solid) line shows our (VMB) result in (l52l . while the blue (dashed) line shows 
the OK prediction (l53l . Clearly the two results agree with each other, as well as with the simulations, in the Gaussian regime 
1/2 < x < 3/4. However, in the outer region 3/4 < x < 1, while the VMB result in (l52l is in perfect agreement with the 
simulation results, the OK result in d53l deviates widely from them. This proves conclusively that in the regime 3/4 < x < 1, 
the VMB result (152t is correct and the OK result (153t is incorrect. 

There is another way to see that the OK asymptotics can not be correct. Since OK theory stems from the asymptotic analysis 
of the following integral representation for the probability distribution of conductance (see eq. 22 of IB9IP : 

yr\G, N) = ™£Ji r**-"^ ( 5 6) 

r(l/8)VW A 7 /syTT2A 
where rj = 2(G/N) — 1, we can now compute from (TS6b the same observable we considered here, namely: 

logF N (Np;G) 

Jg(P) = ~ Jim (57) 

for p 6 K and p ~ 0(1) for N — > oo (see equations ( fTTI i and (fT8T>). and thus compare once again the large N predictions of 
VMB and OK theories against i) numerical simulations, and ii) the exact finite N result, which is fortunately available (see fig. 

Inserting (156b into the definition of the Laplace transform (fTTT i. we obtain after simple algebraic steps: 



2 ( P P 2 



ATI/4 

oo e -Af 2 A(l-p 2 /16) 



x / d\ w& (58) 



n 



<S>n{p) 

The integral 4>n(p) clearly does not converge for \p\ > 4, with the consequence that OK integral representation d56i l fails to 
reproduce the tails of both: 

• Montecarlo simulations of Jq (p) (see fig. [6] left panel), 

• the exact finite N result |[39ll for the Laplace transform in terms of a Hankel determinant (see fig. [6] right panel): 

/v - r 1 ro- + 2)r 2 (. ? + i) 
<w : i 

z=Np 



F N (Np;G) = — det 

c N 



(-d z y +k - 



rd + N + i) m 



which evidently do exist and are perfectly captured instead by our approach in both cases. 
Within the range of validity \p\ < 4, the integral 4>n{p) can be evaluated and gives eventually: 

exp(-N 2 (p/2-p 2 /32)) 
F N (Np;G)= Pi (1 _J /16)1 V » , ifb|<4 (60) 

Note that all TV-dependence has completely dropped out from the prefactor, leaving us with two perfectly compatible (even 
though apparently discordant at first glance) consequences: 



1 . According to OK theory, the limit 



MP) = Hm l0gF ^f P;G) = P/2 p 2 /32 for |p| < 4 (61) 
which is correct but incomplete (as their integral has nothing to say about the tails \p\ > 4). 
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2. Not surprisingly then, from the OK theory, the leading (1/N) order of the cumulants of the distribution, computed through 
the formula (1211 : 



1 



I: 



N J dp' 



logF N (Np; G) 



p=0 



(62) 



is correctly reproduced (see also [38] for an independent calculation of such cumulants, which is in perfect agreement with 
OK result). 

In summary, the OK integral representation in ( TSoT ) is only adequate around the Gaussian peak (the p = neighborhood in 
Laplace space, which is exactly the only region which cumulants probe (see (l62l)). and in this neighborhood it has the merit 
of producing the correct leading 1/N term of the expansion of such cumulants (unattainable by our method) and confirmed 
independently in l38tl . Outside this region, however, the OK integral representation is invalid and the asymptotic analysis of it 
outside its range of validity obviously produces an incorrect result. Note also that, for any fixed N (however large), one has from 
(l56j l that yff^fG = 0, N) = ?c? K) (G = N. N) > strictly, while it is well-known that the density must vanish identically at 
the edges Q7i 13a, l52ll for any N. This simple observation rules out the claims of exactness of d56i l in ll39tl . 



E. Final results for the conductance case 
To summarize, the density of eigenvalues (solution of the saddle point equation (l33l ) has the following form: 



2^T(l-T) 



4+p _ rri 

2p 



2-rrVT V P 



T 



2tT\ 



i=yr-(i-4/b|) 



< T < 1 -4<p<4 

< T < 4/p p>4 
l-4/|p|<T<l p<-4 



(63) 



One may easily check that g* (T) is continuous at p = ±4, but develops two phase transitions characterized by different supports. 
The action at the saddle point is given by: 



-§2+f + 21og2 -4<p<4 



S[g* p } = \ 3/2 + logp p>4 
3/2+p + log(-p) p<-4 



(64) 



which is again continuous at p = ±4. 

Using the above expressions for the saddle point action and the result in d45l l. the expression for the free energy difference 
J G (p) = S[g*} - S[q$\ follows from ©: 



p i_ p 

32 ~ 2 



-4 < p < 4 

Jg(p) = { 3/2 + log(p/4) p>4 
3/2 + p + Iog(-p/4) p<-4 



(65) 



Using this expression for J a (p) in the Legendre transform (I29l > and maximizing gives the exact expression for the rate function 



*g{x) = { \ - \og{Ax) 



1 < T < 3 
4 - X - 4 

< X < T 



i - log[4(l - x)\ l<x<l 
From this formula, one can derive the leading behavior of the tails of Tg(G, N) as: 



s-> in 

J'g(G, N) w exp 



? G (G,N) G 3 W exp 



^N 2 [- log (4G/N)]\ =G? N2 / 2 
N 2 [- log (4(1 -G/N))} 



K T 2 



(N - Gf N2 ' 2 



(66) 



(67) 



(68) 
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in agreement with ||37L|38|,l52ll. 

The most interesting feature of ( f66b is the appearance of discontinuities in higher-order derivatives at the critical points: more 
precisely, the third derivative of ^g^) is discontinuous at x = 1/4 and x = 3/4. In fact, Sommers et al [37] found that for 
finite iVi, N2 there are several non-analytical points. Only two of them survive to the leading order N — > 00 and, in our picture, 
these correspond to a physical phase transition in Laplace space. 

In summary, the following exact limit holds: 



lim 

AT-S-OO 



2 logy g{Nx,N) 
(3N 2 



* G (x) 



(69) 



where the rate function ^q{x) is given in 
The free energy difference in Laplace space 



for large N has been compared with: 



1. Montecarlo simulations over the range p — (—30, 30), which already for N = 5 show an excellent agreement (see fig. [6] left 
panel). For a given p between —30 and 30, the numerical Jg(p) (and analogously for Jq(p)) is computed as: 



Jg(p) 



\ 



rij<fc \ x j -x k \ 2 ) 



(70) 



N 



where the average is taken over N random variables {xi} drawn from a uniform distribution over (0, 1). 

, the exact finite N result from l39ll for the Laplace transform of the density in terms of a Hankel determinant (see d59l). In fig. 
H](right panel), we plot our asymptotic result Jg(p) (l65l l together with — logF jv(-/Vp; G)/N 2 , where Fjy(Np; G) is the exact 
finite N result ((59]l for the Laplace transform from |33l, for -30 < p < 30 and N = 4, 10. Already for N — 4, our Jq(p) 
reproduces the exact finite N formula with an accuracy of two decimal digits. 



V. DISTRIBUTION OF THE SHOT NOISE 

The dimensionless shot noise is defined as P — Yld=i ( 1 — ') Htl • It is convenient to rewrite it in the form P = 
N/4 - E^IiC 1 / 2 - T *) 2 = N / A - Q- The probability distributions of P and Q are related by: 

7 p (P,N) = 7 q (^-P,n) (71) 



^ 4 

It is also necessary to make the change of variable in the joint pdf (O /i, = 1/2 — Tj, so that — 1/2 < fa < 1 /2. The joint pdf 




FIG. 8: Density of the auxiliary Q (schematic). 



16 



(0 expressed in terms of the new variables fa reads: 



N 



p(jj, u . . . , fj, N ) = a n yi \fij - MA/n 

j<k i=l 

and we are interested in the large N decay of the logarithm of J'q(Q, N), where Q = J^i Mi ■ We have: 



V Q (Q,N)=A N 



Again, taking the Laplace transform and converting multiple integrals to functional integrals we obtain: 

poo P 

J Vq{Q-, N)e~^ NpQ dQ — An J D^e"^^ 
where for notational simplicity we keep the same symbols g p and S as before. Of course, the new action S reads: 



(72) 



/ 1/2 " ■ / 1/2 dfii ■ ■ ■ dtXN exp f f X log \& ~ + (f ~ : ) X log Q ~ ^) j 5 (x ~ 



(73) 



(74) 



,1/2 ,1/2 ,1/2 I" ,1/2 

S[q p ]=P g p (fi)fi 2 dfx- / dfidfx'g p (fi)Q p (fj,')log\fj,- fx'\ + C / g p (n)dn-l 

J -1/2 J — 1/2 J — 1/2 J -1/2 

where C is the new Lagrange multiplier enforcing the normalization of the charge density to unity. 
The stationary point of the action S is determined by: 



5g p 







yielding: 



,1/2 

PV 2 + C = 2 d/x'^C/xOlogl/x-M'l 

J-l/2 

Taking one more derivative with respect to /i, we get to the following Tricomi equation: 

P f 1/2 eW) , , 

7-1/2 M-M 

In terms of the solution of ( l78l l. the action (T75T > can be simplified as: 



(75) 



(76) 



(77) 



(78) 



(79) 



where the value of the constant C is determined from ( TTTb by attributing a value to fi within the support of the solution. As in 
the conductance case, we can write the asymptotic decay of Q as: 



y Q (Q,N)e-? N rQdQ: 



exp 



N*{S[g* p ] - S[g* ]} 



exp 



N 2 Jq(p) 



(80) 



Again, in order to solve d78l ) we need first to foresee the structure of the allowed support for f?p(/i). This time, the symmetry 
constraint g*((j) = g*(—fa) reduces the possible behaviors of £>*(/x) to the following three cases: I) g*(fj.) has compact support 
[-L, L] with L < 1/2, or II) g p {n) has non-compact support (—1/2, 1/2), or III) the support of g p (fi) is the union of two 
disjoint semi-compact intervals (—1/2,— L]\J [L, 1/2) with L > (see Fig. [8}. We analyze the three cases separately. 
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FIG. 9: (Color online). Density of shifted transmission eigenvalues fi for N = 6 and p — 12 (theory vs. numerics) for the shot noise case. 



A. Support on [ -L, L] with L < 1/2 
The general solution of (F78b in this case is given by: 

g ;0»)= JL_ =[ C1 -// 2 ] (8i) 

The constant ci is clearly determined as ci = L 2 by the condition that p*(±L) = 0. Thus, the solution within the bounds 
[-L, L] with L < 1/2 is the semicircle: 

p» = ^/L 2- ^ (82) 
where the edge point L is determined by the normalization condition g*(/j,)dfi = 1. This gives: 

L = J^, L <l/2^p>8 (83) 

So, eventually (see Fig. |9): 



w=fv;-^ -vl- M -v? p>8 (84) 

Evaluating the action d79l (after determining the constant C from (1771) by putting /j, = there) we get for p > 8: 



%iS] = j + 5lo g 2 + ilogp (85) 



From @, the value of J Q (p) = S[g*] - S[q$\ (still using S[q%] = 2 log 2) for p > 8 is given by: 

J QW = i + 5 lo g(|) (86) 
Again, the rate function 5'q(.t) is given by the inverse Legendre transform of d86*b . i.e.: 



*q (x) = max [-xp+J Q (p)] = - - 2 log 2 - - log x (87) 
P 4 z 

valid for < x < 1/16. From (ITU , we have the following relation among the rate functions for Q and P: 

*pW = *«Q-«) (88) 

implying for $p(s) the following expression: 

g/ F (a;) = i-21og2-ilogQ-^ for ^<^<J (89) 

In fig. [9] we plot the theoretical density of shifted transmission eigenvalues together with Montecarlo simulations for N = 6 
and p = 12. 
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B. Support on (-1/2, 1/2) 
The general solution of (F78b in this case is given by: 



TTv/1/4-^ 



(90) 



The constant b\ is determined by the normalization condition J_^, 2 0p(^)d^ = 1- This gives b\ = l/p + 1/8. In turn, the 
positivity constraint for the density implies —8 < p < 8. We then get: 



Evaluating the action (1791 ) gives for — 8 < p < 



1 1 2 

P 8 



-2<^<2' -8<P< 



From ©, the value of J Q (p) = S[g*] - S[q$\ for -8 < p < 8 is given by: 

J 0(P) = -256 + g 

Again, the rate function x I'q(.t) is given by the inverse Legendre transform of d93K i.e. 



^q(x) = max \—xp + Jq{p)] — 64 I x - 



8 

valid for 1/16 < a; < 3/16. From the relation d88l ). we have for , J , p(a;) the following expression: 

2 



* P (x) = 64 - - x 



1 



f 1 <r <r 3 
fOT 16- X -16 



(91) 



(92) 



(93) 



(94) 



(95) 



In fig. [T0]we plot the theoretical density of shifted transmission eigenvalues together with Montecarlo simulations for N = 6 
and p = 1. 
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FIG. 10: (Color online). Density of shifted transmission eigenvalues /i for N — 6 and p = 1 (theory vs. numerics) for the shot noise case. 



C. Support on (-1/2, -L] U [L, 1/2) 

For large negative values of p, we envisage a form for the charge density as in fig. [8](rightmost panel), i.e. on a disconnected 
support with two connected and symmetric components. This is because for large negative p, the external potential ppi 2 in (|77| > 
tends to push the charges to the two extreme edges of the box and 1, creating an empty space in the middle. Since we expect 
to have a disconnected support, we cannot directly use the single support Tricomi solution OTb - We need to proceed differently. 
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FIG. 11: (Color online). Density of shifted transmission eigenvalues fi for N = 4 andp = —9 (theory vs. numerics) for the shot noise case. 



We start by recasting eq. d78l ) in the following form: 



1/2 M-M' 



7 d^ + Pr 



pfi = Pr 



1/2 M - M 



1/2 ^(m'K , 
l M - M 

1/2 ^(m'K , 
— ^ 



/x > 
fi < 



In the rhs of d96l ) (first integral) we make the change of variables // — > — fjf, getting: 



pV P /d//+Pr 



Exploiting the symmetry Q*(fx) — g*(—fx), we get: 



,1/2 

= Pr J d(J,'g*(fi') 



1/2 Q*M), , 
~d[i 



(96) 
(97) 

(98) 



1 



1 



M + M' A* — M' 



2 „P t 



2 /? 



(99) 



Making a further change of variables fi 2 = y, [i 12 = y' we get eventually a Tricomi equation for g p (n) = gp(y/ji) / y/JI as: 



i? y-y 



p = Pr 

Solving ( 1 1001 by the standard one support solution (f3Tb and converting back to g*(fi) we get: 

|jPMl(o2 - fJ 2 ) 



7rV(l/4-M 2 )(M 2 -i 2 ) 



(100) 



(101) 



where 02 is an arbitrary constant, fixed by the condition g*(±L) — (the density is vanishing at the edge points). This gives 
a2 = L 2 . Imposing the normalization condition, we get L 2 — 1/4— 2/\p\. The condition that L > implies that this solution 
is valid when p < — 8. Thus for p < —8, we then get: 



MvV-i/4 + 27h 



(102) 



Note that, when p — >• — 8 from below, the equilibrium solution ( 1102t smoothly matches the solution ( f91~b in the intermediate 
regime. The action is readily evaluated from d79l ) as: 



S[Q* P ]=P I d^ 2 g*( f i)~2 dfig* p {[i) log fi 

J L JL 

= j + ilog(b|)-M + ilog2 



(103) 
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The corresponding Jq (p) is given by: 



Jq(p) 



!i f\p\ 
2 log 



M 

4 



(104) 



from which the rate function \&q(x) can be easily derived: 



* Q (V) = max[-xp + J Q (p)} = j -21og2- ^ log f- - x 



3/16 < a; < 1/4 



(105) 



In fig. QT|we plot the theoretical density of shifted transmission eigenvalues together with Montecarlo simulations for N = 4 
and p — —9. 

D. Final results for the shot noise case 
To summarize, the density of the shifted eigenvalues {^} (solution of the saddle point equation d78l ) has the following form: 



f?p(^) 



ir y P 
p_ 



< M < 



P > 8 



8+£ 

Sp 



— 1/2 < /J. < 1/2 -8<p< 



(106) 



bMlVM 2 -l/4+2/|p| 



-1/2 < M < - VV* - 2 /W v VV4 - 2/|p| < M < 1/2 P < -i 



One may easily check that g*(/i) is continuous at p = ±8, but develops two phase transitions characterized by different supports. 
The saddle-point action (T79b is given by: 



( 3 



%*] = < 



i + ilog2 + |logp 
|-^ + 21og2 

M 



p > 8 

-8 < p < 8 



(107) 



[|+ ilog(|p|)-f + |log2 P <-8 



which is again continuous at p = ±8. 

From J28l ), the expression of Jq(p) = S[g*) — S [qq] is: 



Mp) 




V > 8 
-8 <p< 

|log(^)-^ P<-8 



(108) 



from which one can derive (in complete analogy with the conductance case) the rate function for the auxiliary quantity Q: 



i-21og2- ilogz 



®q(x) = <64 (x-IY 



< x < 1/16 
1/16 < x < 3/16 



(109) 



|-21og2-|log(i- : r) 3/16 < x < 1/4 
and from the relation $p(a;) = ^q^I/A — x) one readily obtains the rate function for the shot noise in ( TTOb . 



VI. DISTRIBUTION OF MOMENTS T„ FOR INTEGER n 

In this section, we deal with the more general case of integer moments 7 n — J2iLi Tp, in particular focussing on the case 
n = 2. The conductance is exactly given by Ti while the shot noise is Ti — T2. While we could use the general method outlined 
in Section Hn] with the choice a(T) = T n , as was done for the conductance case (n = 1), it turns out that one can obtain the 
same final results by using a short-cut which combines, in one step, the saddle point evaluation in ( 1261 ) and the maximization of 
the Legendre transform in (29) . Of course, both methods finally yield the same results, but this shortcut explicitly avoids any 
Laplace inversion. Here, we illustrate the short-cut method for the case a(T) — T n , but it can also be used for other linear 
statistics. 

The distribution of the moments 3V„ (%i — Nt, N) is given by: 
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3V„(T„ =Nt,N) = A 



jf 1 • • • jf 1 dT x ■ ■ ■ dT N exp U J2 log IT, - T k \ + (| ~ l0 8^j 6 (E T " - N *j 

(110) 

The short-cut consists in replacing the delta function by its integral representation: 5(x) = J ^e px where the integral runs in 
the complex p plane. The rest is as before, namely that in the large N limit, one replaces the multiple integral by a functional 
integral introducing a continuous charge density g p (x). This gives 

9V„ OTn = Nt, N) = A N J ^-D[g} e -l N2s M (111) 

where the action is given by: 



S[Qp\ = P ( / Q p {x)x n dx — t I — / / dxdx' q p (x)qp(x') log \x — x'| + C I g p (x)dx — 1 

\./o / "'o Uo 



(112) 



where the rhs of d 1 1 II ) is now extremized with respect to both g(x) and p. Notice that here we have already performed the inverse 
Laplace transform of d23l ). Hence the two methods are exactly identical. 
Extremizing the action gives the following saddle point equations: 



px n + C = 2 dx'g*(x') log \x - x'\ 



dxx n g*(x) 



(113) 
(114) 



which in turn determine pas a function of t. 

Multiplying ( II 13b by g*(x) and integrating over x, the action at the saddle point can be rewritten in the more compact form: 



S[q*] 



[ P t + C] 



(115) 



where, as before, the constant C has to be determined from dl 131 ) by using a suitable value of x which is included in the support 
of the solution. For large N, ( 111 II ) gives 



3V n C« = Nt, N) s» exp 
= exp 



4n 2 {S[q* p }-SIq* }} 



where the rate function is given by 



*t„(<) = - s[eg] = s[ e * p ] - 2 log 2 



(116) 



(117) 



having used again S[qq) — 2 log 2 from d45] >. 

Upon differentiation of ( II 131 ). we obtain the Tricomi equation for g p {x): 



n^-x^ 1 = Pr 
2 



g*Jx')dx' 



(118) 



to be solved for different supports of g p (x) as one varies the argument t and consequently the parameter p. As usual, depending 
on the value of p, we need to first anticipate the 'type' of the solution, i.e., the form of its support and then verify it a posteriori, 
as illustrated below. 
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A. Large p: support on (0, L p 



Consider first the case when p is very large. Since the external potential in (II 131 ) is of the form px n + C which is rather 
steep for large p, we anticipate that the charge fluid will be pushed towards the left hard edge at x — 0. In this case, the general 
solution of dl 181 ) with the restriction n £ N is given by: 



1 



ir 2 ^fx(L p - x) 

Lp y/x'(L p - x^npix') 71 - 1 dx' 



2 



x — x 



C x 



where C\ is an arbitrary constant. Evaluating the principal value integral yields 

1 



Q* v (x) = - 



Tr 2 y/x(L p - x) 
npL"y/WT{n- 1/2) 



F\ l,-n;3/2-n;— + d 



4r(n+ 1) 

where 2-^1(0, b; c; x) is a hypergeometric function, defined by the series: 

, , ab a(a+l)b(b+l)x 2 

2 F X (a, b- c; x = 1 + -x + [ > \ ' — 
c c(c + 1) 2! 

Determining the constant C\ by the requirement that g*{L p ) = 0, we obtain: 

1 



g P (x) = 



7r(2n — l)y/x(L p — x) 
x [ 2^1(1, -n; 3/2 - n; x/L p ) + 2n - 1] 



(119) 



(120) 



(121) 



(122) 



The edge point L p is finally determined by the normalization requirement J a p g*(x)dx = 1, yielding after some elementary 
algebra: 



L p = (2V5?r(n)/(pr(n + 1/2))) 



l/l! 



(123) 



Imposing now ( II 14b as: 



7r(2n - 1) 



y/x(L p - x) 



[2n - 1 + 2^1 (1, -n; 3/2 - n; x/L p )] 



(124) 



we obtain a very simple relation between p and t: 



tn 



(125) 



Armed with (1125b and ( 11221 ). we can now evaluate the action ( II 15b eliminating p as: 



= -L + 21og2 + -log 
In n 



r(n+l/2) • 
2^r(n + 1) 



(126) 



Equation (11221 ) is valid as long as L p < 1 (the edge point of the support such that g p (L p ) = 0). From ( 11231 ), putting L p = 1, 
one thus finds that the solution is valid for p > p\ where 



Pi 



r(n+l/2)' 



(127) 
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Consequently, from ( 1125b . it follows that the solution is valid for t < t*, where 

r(n + l/2) 



2V5FT(n + l) 



(128) 



As a check, for n = 1 we have t* = 1/4 from (1128b . L p = A/p = At from (1123b and (1125b and the condition L p < 1 implies 
i < t* = 1/4 as expected (compare with subsection lTV Ab . Forn = 2, we have t* = 3/16 and p* = 8/3. Forn = 2, this regime 
(t < t* = 3/16 and hence p > p* — 8/3) corresponds to the leftmost panel in Fig. [16] 

In this region of t, the rate function is easily computed as 



*T n («)=%*]-21og2 
1 



1 / T(n+l/2) 

2 + ° g ^2ViT(n + l) 



-lost 



(129) 



Combining ( 1129b with ( 11 16b . one obtains as a new result the precise left tail asymptotics for the n-th integer moment distribution: 



y Tn (7 n = Nt,N)nexp 



/3N 2 
In 



log i =i ! » 



(130) 



In fig. Q~2]we plot the theoretical density of eigenvalues for the n = 2 case together with Montecarlo simulations with A = 5 
and p = 5. 




FIG. 12: (Color online). Density of eigenvalues x for TV" = 5 and p = 5 (theory vs. numerics) for n = 2. 



B. Intermediate p: support on (0, 1) 



We now look for the solution of 



rZx n - x = Pr 



■dx' 



(131) 



with a nonzero support over the full allowed range [0,1]. Using the general solution in ( T3"T1 ) with the choice a = and b = 1 we 
get 



npr(n- 1/2) 
AT(n + 1)tt 3 / 2 



2 Fi(l,-n;3/2-n;x)-C„ 



^(1 - at) 



(132) 



where C n is an arbitrary constant to be fixed by g*(x)dx = 1, yielding 



4y^r(n) 

pr(n- 1/2) 



(133) 



so eventually: 



pT(n- 1/2) 
4r(n)7r 3 / 2 



a Ji(l,-n;3/2-n;s) + 



4y?r(rt) 

pI>-l/2) 



(134) 
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Next, we need to impose the condition ( II 14b : 



t = 



dx x n g*{x) 



leading to: 



t 



T(n+ 1/2) 



1 



pT(n - 1/2) 



v^r(n + 1) 

From ( I136l l we can derive the relation between p and t as: 

p = a n - b n t 

where: 

4^r(n) 
fln " r(n + l/2) 
4?rr(n)r(n- 



(1 - 2n) 



&n = 



1 



[r(n + l/2)] 2 

Inserting ( 1137b and (1134b into dl 15b . we obtain, after a few steps of algebra, the action: 

* Le p J = — t - 



r(n + 1/2) 



0Fr(n + 1) 



log 4 



(135) 

(136) 
(137) 

(138) 
(139) 

(140) 



Next we need to determine the range of validity of this solution. This is obtained simply by the fact that the density g* (x) in 
132b must be positive. Let us first rewrite the solution (1132b as 



Q*Jx) = 



1 



A p (x) 



where 



P r(n-l/2) 

^ (x) = 1+ 4^r(.) 2 ^ l(1 '' 



; 3/2 — n; x). 



(141) 



(142) 



To ensure £>*(:e) > 0, we have to just ensure that A p (x) > in (11421) . How does A p (x) vary as a function of a: in a; G [0, 1]? 
It can be easily seen that this function has a global minimum at some intermediate value < x* < 1. To ensure its positivity, 
we then have to ensure that A p (x*) > 0. This will be true only for a range of values of p, i.e., when p\ < p < p* (where p\ 
is precisely the lower edge of the validity of regime I in the previous subsection and is given in ( 11271 )). Consequently, using 
( |137b , this sets a t range t* < t < i£ f° r me validity of this regime II, where t* is given in ( 1128b . Now, for arbitrary n, p\ and 
consequently t\ have rather complicated expressions which we do not detail here. But for n = 2, their expressions are rather 
simple and we get 



16 



and <2 



(143) 



This then defines regime II with a full support over [0, 1], namely —16/3 < p < 8/3 and consequently 3/16 < t < 3/4, is 
shown as the second (from the left) region in fig. [16] 

Thus in this regime II where t* < t\, the rate function, for arbitrary n, has a quadratic form: 



thus implying (fT3l : 



*T n (t) = S[Q* p ] 



T T „(T„,^)«exp 



1 A ^ 

log4 = — 



2 2 



T 

j n 

~N 



T(n+ 1/2) 
V^T(n + 1) 



r(n + l/2) 



1 2 



0FT(n + 1) 



(144) 



(145) 
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FIG. 13: (Color online). Density of eigenvalues x for N — 6 and p = 1 (theory vs. numerics) for n = 2. 



In fig. Q~3]we plot the theoretical density of eigenvalues for the n = 2 case together with Montecarlo simulations with TV = 6 
and p = 1. 

From the Gaussian shape in ( 11451 ). the mean and variance of T n can be read off very easily: 



NT(n+ 1/2) 

V^r(n + i) 



var(T„) = — 



[r(n + l/2)] 2 



2^7rr(n)r(n + 1) 

Novaes j25ll recently computed the average of integer moments and obtained the following expression for arbitray n £ N 

var (T n ) 

. 08f 



(146) 
(147) 
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FIG. 14: (Color online). var(T„ ) as a function of n G N dT4Tb for /3 = 2. In red the asymptotic value «* = 1 /2vr/3 w 0.07957 ( TT49t . 



■CJ„> =: V(^]4 " 



(148) 



Using elementary properties of Gamma functions, it is easy to show that formulae ( 11461 ) and (11481 ) do indeed coincide, a fact not 
completely apparent at first sight. 

Conversely, the exact expression for the large N variance ( 1147b is new 4 , since the general integral in Beenakker's formula 
l29tl does not appear easy to carry out explicitly. Obviously, d!47t agrees with the known result [29] for the conductance 
var(G) = l/8£ for n = l. 

From ( 11471 ). it is easy to extract the asymptotic value: 



v* = lim var(T„) 
n— >oo 

which is plotted in fig. [l4]for f} = 2 together with ( 1147b . 



1 

2/3^ 



(149) 



Novaes has recently computed the variance for any finite number of open channels Ni and 7V2 f 25] : however, extracting the asymptotics from his formula 
does not appear to be easy. 
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For general n, as one increases the value of t, so far we have seen two regimes: regime I (0 < t < tj) with support over [0, L p ] 
and then regime II (t* < t < i£) with support over the full range [0, 1]. What happens when t increases beyond t^l For arbitrary 
n, the analysis becomes rather cumbersome. So from now we restrict ourselves only to the case n = 2 (which turns out already 
to be rather nontrivial). But at least for n = 2 we are able to obtain a full picture and in the two subsections below we show that 
apart from regime I and regime II already discussed above, two further regimes appear as one increases t beyond t$: regime III 
(for t\ < t < £3) where the solution has a disconnected support with two connected components discussed in subsection I VI Dl 
and regime IV (for t$ < t < 1) where the solution again has a single support but on the other side of the box over [M p , 1]. Since 
the solution in regime IV is simpler (single support), we will first discuss this case in the next subsection I VI CI and finally the 
more involved case of regime III (with a disconnected support) will be discussed in subsection lVIDI 



C. Support on [M p , 1) (n = 2): regime IV 

Focussing on the n = 2 case, we now look for a solution of dl 18b with a single support [M p , 1] where M p is yet to be 
determined. Using the general single support Tricomi solution (l3TT l choosing a = M p and b = 1, one obtains the following 
explicit solution 



Qp( x ) 



. 1 J(l-x')(x' - M p ) , , 
p Pr / ^ ^ ^ x' dx' + D 

Mr, X — X 



n 2 y/(x-M p )(l-x) 

where D is an arbitrary constant. Evaluating the principal value integral in d 1 50b and imposing g*(M p ) = we obtain 



(150) 



e>) = -^ + M,- 1} ^ (151) 

The lower edge M p is determined by the normalization condition J M g*(x)dx — 1, yielding a quadratic equation for M p : 

p(M p - 1)(1 + 3A/p) = 8 with two roots M p = (1 ± 2^1 + 6/p)/3. Noting that when p -00, it follows from physical 
consideration that the charge density must be pushed to its rightmost limit indicating that M p — > 1 as p — > —00. This condition 
forces us to choose the correct root as 



The condition < M p < 1 implies for p the condition p < p\ = — 6. The relation between p and t is then obtained using dl 141 ), 
resulting in the condition: 

15 + 27M„ + 13M? + 9M^ 

7 =t (153) 

16(3M P + 1) 

where M p is expressed as a function of p in (11521 ). The solution of ( 1 153b is quite cumbersome to write down explicitly, but is in 
principle feasible. Note that when p — > —6 from below, M p — >■ 1/3 from (11521 ) and consequently from d 1531 ). t — > £3 = 29/36 
from above. In other words, the solution (1 1 5 1 b is valid in regime IV defined by 

29 

P < Pa, = —6; consequently, t% = — < t < 1 (154) 

36 

This regime IV is shown in the extreme right part of Fig. [16] 

Once p has been determined has a function of t from (11521 ) and ( 1 153b and substituted into the density d!5U , the action and 
the rate function can be computed from d 1 15b and dl 17b , by evaluating numerically the corresponding integrals. We omit these 
details here. 

In fig. [I5]we plot the theoretical density of eigenvalues for the n = 2 case together with Montecarlo simulations with TV = 4 
and p = —10. 



D. Disconnected support (n = 2): regime III 

For n = 2, it then remains to find the solution of dl 18b in the narrow band f% = 3/4 < t < t% = 29/36 or equivalently for 
£>3 = —6 < p < P2 — —16/3 (see Fig. [16}. This is the regime III. Let us first try to anticipate what the solution may look like in 
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FIG. 15: (Color online). Density of eigenvalues x for N — 4 and p — — 10 (theory vs. numerics) for n — 2. 



this regime. For this, let us consider the two regimes namely regime II and regime IV respectively to the left and right of regime 
III. 

Consider first regime II (t* — 3/16 < t < t\ — 3/4). In this regime the solution g*(x) has a single support over the full 
range x £ [0, 1] given in ( 1134b . which for n — 2 (using the special value of the hypergeometric function) simply reads, with 

P* 2 = -16/3 <P<Pl = 8/3, 



1 



TTy/x(l — x) 



P 



l + |(l + 4x-8^) 



Now, when t tends to the maximum allowed value in regime II, namely, t — s- t\ 
— 16/3 from above, the solution in ( 1155b tends to 



(155) 



3 /4 from below or equivalently p —> p\ 



16/3 



(x) 



16 



3lTy/ x(l — X) 



(156) 



with a quadratic minimum at a; = 1/4 where the density vanishes. This is just the edge of regime II. If t increases slightly 
beyond t% = 3/4, this single support solution is no longer valid. However, it gives the hint that for t > i| — 3/4, the charge 
density must separate into two disjoint components, one on the left side over [0, l\] and one on the right side over [Zq, 1] with an 
empty stretch [Zi , I2] separating them. This empty stretch must increase as one increases t beyond t% = 3/4 in this regime III. 
Indeed, as t increases further, the left support [0, Z1] must shrink in size and the right support [Z2, 1] must increase in size and 
finally when t hits the value = 29/36, l\ must shrink to and Z2 must approach M p = 1/3 and then one arrives in regime IV 
discussed in the previous subsection. At exactly t = t$ — 29/36 or equivalently at p = = — 6, the solution (border of regime 
W) can be read off (ED with M p = 1/3 



7T % /T _ ¥ V 3/ 



.3/2 



(157) 



Thus, the solution in regime III, namely for p\ — —6 < p < p\ — —16/3 must interpolate between the solutions given 
in J 156b and in ( 1157b valid respectively at the two edges of regime III and have a disconnected support with two connected 
components over [0, Zi] and [Z2, 1]. We were able to find this solution explicitly. Its derivation is outlined in the Appendix. The 
result reads 



8p(x) 



-P 



y/(x-h)(x-l2) 3 



ir^/x(l — x) 

which is valid for all < x < l\ and I2 < x < 1 and the two edges Zi and Z2 are given by 



(158) 



h 

h 



l-3i/l 



16 



16 
3^ 



(159) 
(160) 



Note that Zi and I2 are real only if p < — 16 /3. Furthermore 1% > only if p > —6. Thus this solution is valid over the full range 
j>3 = —6 < p < P2 = —16/3. This then defines regime III. Note also that the solution ( |158b smoothly interpolates between the 
solutions in (1 156b (whenp — > —16/3) and in d 1 57b (whenp — > —6). 
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The relation between t and p can be obtained substituting ( 11581 ) in ( II 14b and performing the integral. Similarly the action and 
the rate function can be computed by using the exact density and performing the integrals in dl 151 ) and dl 171 ) by Mathematica, the 
details of which we omit here. The Montecarlo simulations in this regimes are much harder to obtain due to large fluctuations 
in sampling from the leftmost residual band < x < l\ and a very large N is necessary to achieve a satisfactory picture. 
Nevertheless, we observed upon increasing N a trend in the equilibrium density which is fully compatible with the analytical 
disconnected-support solution found above. 



E. Final results for the n = 2 moment case 

In fig. [16] we propose a schematic summary of the different phases of the density of integer moments for the case n = 2. 
As a consequence, the rate function will display three (and not just two, as for conductance and shot noise) non-analytical 
points corresponding to physical phase transitions in Laplace space. Starting from high p values, the fluid particles are initially 
confined towards the left hard edge. Then, when p hits the critical value p = 8/3, the fluid spreads over the entire support [0, 1]. 
In the narrow region —6 < p < —16/3 the density splits over two disjoint (non-symmetric) components of the support, and the 
leftmost disappears upon hitting the value p — —6, leaving the charges leaning against the right hard wall. 




FIG. 16: Schematic table summarizing the relevant phases in the t and p space for the density of n — 2 moment. As t varies in the allowed 
range [0, 1] and consequently the parameter p in [— oo, oo], the fluid density displays 4 different phases as shown in the top panel. These 4 
regimes are separated by three critical points: tl = 3/16 (consequently pi — 8/3), t% — 3/4 (p£ = —16/3) and t%, — 29/36 (P3 = —6). 
Consequently, the rate function ^fT n (t) has 4 different expressions according to different regions in the t 6 [0, 1] segment. 



VII. CONCLUSIONS 



The Coulomb gas analogy, together with recently introduced functional methods, brings about an efficient formalism for the 
computation of full probability distributions of observables (valid for a large number of electronic channels in the two leads) in 
the quantum conductance problem. Generically, the distribution of any linear statistics of the form A = YliLi a (Ti), where the 
Tj's are transmission eigenvalues of the cavity and a(x) is any smooth function, can be derived within the formalism described 
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in this paper: the problem amounts to finding the equilibrium configuration of an associated 2d charged fluid confined to the 
segment [0, 1] and subject to two competing interactions: the logarithmic repulsion, generated by the Vandermonde term in the 
jpd (0, and a confining potential whose strength is tuned by the Laplace parameter p. Interestingly, this auxiliary Coulomb gas 
undergoes different phase transitions as the Laplace parameter is varied continuously, and this physical picture is mirrored in 
the appearance of very weak singularities in the rate functions of observables at the critical points. Already the leading TV-term 
of the free energy of such a gas (the spherical contribution) displays non-Gaussian features in the tails, while the central region 
obeys the Gaussian law, in agreement with the general Politzer's argument ll59ll : conversely the tails follow a power-law decay 
and the junctions of the two (or more, as in the case of integer moments) regimes are continuous but non-analytical points. Note 
that it is not necessary to develop a (1/N) expansion of the free energy and look for higher genus terms to appreciate deviations 
from the Gaussian law. From the central Gaussian region, one can easily reads off the mean and variance of any linear statistics 
of interest: this way well known results for e.g. conductance and shot noise are recovered and new formulas (such as the variance 
of integer moments and its universal asymptotics) can be derived. Our results are well-corroborated by Montecarlo simulations 
both in the real and Laplace space, as well as by comparison with exact finite N results when available, which convincingly 
disprove the large-iV asymptotic analysis performed in J39ll . In summary, the Coulomb gas method (well-suited to large N 
evaluations) reveals a rich thermodynamic behavior for the quantum conductance problem already at the leading order level, 
and we expect that it will enjoy a broad range of applicability. 
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Appendix A: Explicit Two Support Solution for n = 2 in Regime III 



We consider the integral equation (111 81 > for , 



Z 1 Q*Jx')dx' 

px = Pr / pV - (Al) 

o x x 



and look for a solution that has a disconnected support with two connected components ([0, li] and [I2, 1]). This solution is valid 
in regime III discussed in subsection VI-D. While a general single support solution of the singular integral equation can be found 
using Tricomi's formula given in (f3Tb with g(x) — px and discussed in subsections VI-A, VI-B and and VI-C, it is much more 
complicated to obtain an explicit solution with more than one connected component of the support. To find such a solution, 
we actually use an alternative method originally used by Brezin et. al. to find a single support solution of the singular integral 
equation in the context of counting of planar diagrams 16011 . This method consists in making a judicious guess for the solution 
and then uses the uniqueness properties of analytic functions in a complex plane to prove that the guess is right. Although, for the 
single support solution one does not have to use this route since the explicit general solution of Tricomi is available ( the authors 
of ref. 116011 were perhaps unaware of the general single support explicit solution of Tricomi). Nevertheless, this alternative 
method of i60ll can be fruitfully adapted to find a two support solution (as in our case) in a simpler way (as demonstrated 
below), where a general solution is somewhat difficult to obtain explicitly. 

Let us assume that the solution g*(x) of ( IA1 b has a disconnected support with two connected components [0, ii] and [I2, 1] 
where l\ and I2 > h are yet to be determined. Generalizing the route used for single support solution in ref. Il60ll . we first 
introduce an analytic function (without the principal part in ( 1AU ) 



F(x) = [ 
Jo 



1 g*(x')dx' f 1 g*(x')dx' 



(A2) 



defined everywhere in the complex x plane outside the two real intervals [0, li] and [l 2 , 1]. This new function F(x) has the 
following properties: 

1. it is analytic in the complex x plane outside the two cuts [0, l\] and [I2, 1], 

2. it behaves as 1/x when |x| — > 00 since J g*(x')dx' = 1 due to normalization, 

3. it is real for x real outside the two cuts [0, li] and [I2, 1], 

4. as one approaches to any point x on the two cuts [0, l±] and [I2, 1] on the real axis, F(x ± ie) = px =p ing^x). This last 
property follows from dAU . 
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From the general properties of analytic functions in the complex plane it follows that there is a unique function F(x) which 
satisfies all the four properties mentioned above. Thus, if one can make a good guess for F(x) and one verifies that it satisfies all 
the above properties, then this F(x) is unique. Knowing F(x), one can then read off the solution g*(x) from the 4-th property. 
It then rests to make a good guess for F(x). We try the following ansatz for F(x) valid everywhere outside the two cuts [0, h] 
and [Z 2 , 1] 

F(x) = px- P ^(x-l^ix-hf. (A3) 
V x ( x - !) 

This ansatz clearly satisfies the first property. Now, expanding F(x) for large |x| we get 

F(x) -> + + 3l 2 )p+ 

+ |(-3 + 2h + + 6l 2 - 6hh ~ 3I|) - + 0(x- 2 ) (A4) 

Since the second property dictates that F(x) — > l/x, it follows that we must have 

h + 3l 2 = l (A5) 
(-3 + 2h + If + 6h - 6Z1/2 - 3Zl) = 8 (A6) 

Eliminating^ in ( IA6I 1 using iA5\ gives a quadratic equation for U . 2Zf — Zi — 1 = 6 jp with two solutions: l\ = (l±3^/l6/p)/4. 
The correct root is chosen by the fact that when p — > —6, l\ — > as follows from the solution in ( 1157| i. This then uniquely fixes 
li and l 2 given respectively in ( 11591 ) and 1160) . 

The ansatz F(x), with the choices li and l 2 as in J 159b and ( 1160b then satisfies the second property. It is easy to check that 
F(x) satisfies the third property as well. From the fourth property one then reads off the unique solution as given in (11581 ). This 
two support solution is clearly valid only in the regime III namely forp\ = — 6 < p < p 2 = —16/3 and it smoothly matches 
with the solutions of regime II and regime IV respectively asp^ —16/3 and p —> —6. 
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